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ABSTRACT 


A  sea  floor  benchmark  experiment  was  conducted  in  an 
area  about  16  nmi  west  of  Pt.  Lobos,  California  (36°30'N  x 
122°17'  W)  during  18-22  May  1985. 

Two  baseline-crossing  methods  were  used  to  determine  the 
relative  positions  of  acoustic  bottom  transponders.   The 
method  of  least-squares  adjustment  was  used  to  analyze  the 
data.   Relative  position  determination  of  the  transponder 
array    is  discussed  and  recommendations  are    made  for  futher 
improvement.   The  advantage  of  these  methods  is  their 
simplicity.   Their  disadvantage  is  the  relatively  large 
amount  of  ship  time  they  require  to  achieve  acceptable  accu- 
racies. 

Transponder  arrays  such  as  the  one  deployed  can  be  used 
for  solving  many  types  of  problems  in  sea  floor  engineering, 
which  will  be  of  increasing  importance  in  the  future. 
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I-    IH2RQ duct ION 

A-        BACKGROUND 

In   todays   world,    the   ocean  scientist    and  engineer   need 
increased   accuracy   for   marine    positioning    to   accomplish 
their   work.    The    establishment    of   a    system   for   the   precise 
determination    of    positions   on    the   deep   seafloor    is    one   of 
the    most   fundamental    challenges   of   working   in   the  open    sea. 

A   seafloor   benchmark    positioning    system   is    proposed 
which   can   determine    accurately   the   location    of-  objects    on 
the    sea    floor    (Saxeua,    1974).      Such    precisely   located   benca- 
marks   can   later  be  used   in    turn   to    delineate  offshore   prop- 
erty   lines    and   national   boundaries-       They   are    also   useful    in 
solving    problems    in    seafloor   engineering,    in   plate    tectonic 
studies    and  in   connection   with    bore-hole  reentry   associated 
with    off-shore   oil  recovery. 

The    objective   of    this    study  is    to   demonstrate   the   feasi- 
bility  of   establishing  such    seafloor   benchmarks. 

B.       DESCRIPTION 

Early   experiments   associated   witn   seafloor    acoustic 
transponder  arrays   were   carried  out   by    Hart     (1967),    Haehnle 
(19o7),    Fubara   and  Mourad    (1972),    and    McKeovn    (1975),    wno 
used    a   baseline   crossing    method   to   solve    for   the   distances 
between    transponders.      "Baseline"   in    this    tnesis   means    the 
line    between    the    projections   of    two   bottom-mounted    acoustic 
transponders   onto    the   sea   surface. 

The    baseline    accuracies    obtained   by   Hart    (1967),    HcKeown 
(1975),    and    Fubara   and   Mourad    (1972)    were    ±15-7    m   for   base- 
lines  of    5509    m,    ±3. 1    to    ±4.2    m   over    baselines    of    b373    to 
7219    m   and    ±15-5    m   over   baselines    of    9364    m,    respectively. 
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The    first   two    of    the    above   authors   did    not   use    a  least 
squares   method   for   data  analysis. 

There  are   three    parts   to   the   determination    of    the   abso- 
lute   position    of    an    ocean-bottom   transponder   array    using 
satellites-      They   are    (Figure    1.1); 

1.  To  determine    the   geodetic   position   of   tne  ship's 
receiving    antenna   from    the  Global   Positioning    System 
<GPS)  ; 

2.  To   convert   the    geodetic   position    of   the    antenna    to 
the   ship- mounted   acoustic   transducer;    and 

3.  To   determine    the   geodetic   position   of   the   acoustic 
transponder  array   on    the   sea   floor   using    a   snip- 
mounted   acoustic   transducer. 

This   thesis   is  concerned    with   the    relative    position   of    a 
bottom- mounted   transponder   array  with    regard  to    the    ship- 

C.       DESCRIPTION    OF   THE   SYSTEM 

Acoustic    transponders   manufactured   by   Oceano   Instruments 
were    used  for    this   project.       They   were   deployed    by    the    R/V 
Acania.      Figure    1.2    shows   the    configuration   of    the   Oceano 
system   used.    '  Their    principal    features   are: 

1.  AM-121    :     (Acoustic   Module) 

The   hydrophone   was    fixed   amidships   to   the 
port  side    at    a   depth    2. 95   m. 

-  Transmission    frequencies    :    8    to    16    kHz. 

-  Reception   frequencies    :    8    to    16    kHz. 

2.  ST-121    :     (Recoverable    Transponder) 

-  Frequencies    :    8    to    16    kHz. 

-  Delay   time    :    15.0±0.1    ms. 

3.  TT-201  ;  (The  Accoustic  Telecommand  Module) 
has  the  following  functions: 

-  Code  acoustic  signal  transmitted  by  either  the 
acoustic  module  or  UQC  transducer. 

-  Monitors  range  to  desired  transponder- 


1  1 
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Figure    1.  1 


The    Configuration   of    Bottom   Transponder   Array 
Shown   with   Ship   and   GPS    Satellite. 
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-  Monitors   reception    of    code   by    transponder. 

4.  IM-100    :     (Data   Interface   Module) 

-  Accepts    up    to   16   data   inputs. 

-  Max.    converted   data   format    :    7    digits    without 
sign    or   6   digits   plus    sign. 

-  Computer   interrupt    reguest   possible   through 
any   input. 

-  Each    channel   may   be   read   separately. 

5.  RM-201    :    (Rangemeter) 

-  Transponder    turnaround   delay   offset    :    0   to    99.99   ms. 

-  Mean    velocity  of   sound   input    :    0    to   9999.9    units   per 
second. 

-  Standard  reception    filters   ;    8    to    16    kHz    (0.5   kHz    step) 

-  Bandwidth    :     203    Hz. 

Critical  to   the    experiment   are   the    following: 

1.  The  transducer  (a  part  of  the  AM-121)  was  2.95  m 
below  the  surface,  and  the  transmission  and  reception 
frequencies  were  8  to  16  kHz. 
2-  The  system  delay  on  the  RM-201  was  set  at  15  millise- 
conds, i.e.,  the  time  lapse  between  reception  and 
transmission    of  a    signal. 

3.  Velocity    on   the   RM-201    was    set    to    1480    m/s. 

4.  Interrogation    period    was   set    at    10   s. 

5.  Interrogation    frequency    (kHz)     was    set    15    kHz. 

The    numerical   constants    of     (3)    and    (4)     were    used    for   the 
survey.      Details    of    the   transponders   are    given    cy   Oceano 
Instruments    (1984). 
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Configuration  of  The  Oceano  Systems  Used. 
1U 


II-     THE    OPERATIONAL   PROCEDURE    FOR    IMPLANTATION    OF    A 

TRANSPONDER   ARRAY 

A-       SITE   SELECTION 

The    main   criterion  for   the    site   of    the   experiment   was 
that    it    be   relatively   flat    over   an    area    approximately   2    km    x 
2   km.      Since   there  were  no   existing   caarts   detailed    enough, 
a   reconnaissance    survey   was    perfomed   by    the    R/V    ACANIA   on    1 0 
and    11    April    1985.       Some    10    north-south    sounding   lines   and 
1-cross    line   were   run   with    1-nmi   spacings   and  -5-min    position 
fixes. 

The   area   chosen    for   the    experiment    (Figure    2.1)     is 
approximately    24    nmi    west   of   Pt.    Lobos,    California,    in   water 
approximately    1500   to    2000    m    deep.      It    is   the  area   from 
36°29.5«    to    36°31»N    and    122°16«    to    122017.8»W.       Mini-Hanger 
stations    (Figure    2.2)     were    set    up   at   Point   Sur 
(36°18«  20.279"    N    X    121°    53*56.  179"W)     and   Carmel 
(3b°33*  49.  176"    N    X    1 2  1°53 • 48. 358"W)    for    additional    ship 
position  control. 

B.        IMPLANTATION    PROCEDORE 

The    transponders    were   deployed   at    pre-selected    sites 
using   a    buoy-first,    anchor-last,    technique.      The   anchors 
used    were    railroad   wheels    which    weigh    about    290    kg    in   water. 
The    mooring   cables   between    the   anchors   and    the    transponders 
and    between    the    transponders  and   the    buoys   are    each    8-m 
long.      The    transponders   weigh    about    46   lbs   in   water,    and 
each    buoy    has    55    lbs    positive    floatation.       Three    17-in 
Benthos   glass    sphere    floatation    buoys    were    attached    to    each 
transponder    by    means    of    5/32-in   stainless  cables    8-m    long. 
Altnougn    2    spheres    would    have    provided    sufficent    buoyancy 
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for  recovery,  a  third  was  used  to  provide  additional  tension 
(about  50  kg  total)  to  reduce  the  motion  of  the  transpon- 
ders.  The  mooring  system  is  shown  in  Figure  2.3. 


NOT  TO 
SCALE 


3  17-inch  Benthos  glass 
spheres   in  hard  hats  - 
55   lbs  buoyancy  each 


ff 


Oceano  acoustic  transponder 


{'4 


k 


5/32-inch  stainless 
8  m  long 


release  mechanism 


5/32-inch  stainless 
8  m  long 


ail  road  wheel   used  as 
anchor 


OCEAN     BOTTOM 


//////// '  ■ ' / //? 


Figure  2.3    Schematic  of  Bottom  Transponder  Mooring 
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III.  DETEHHINAIION  OF  HARMONIC  MEAN  SOUND  VELOCITY 

A.  I  NX  HOD  UCT  ION 

Historically ,    echo  soundings  have    been   made    by   assuming 
a   constant,    approximately   average,    value   for   the   velocity   of 
sound   throughout    the    water   column,    usually 

4800    ft/s    (=1463    m/s)    or    1500    m/s,    and   then    cocrecting    for 
variations   of    this  assumed    value   for   the    actual    water 
column.      Echo    sounders  are    time-measuring   devices.    The 
sounder's  acoustic   transducer   emits  a    sound    puise    which    is 
reflected   upwards   by    the    bo ttom- mounted    acoustic    transponder 
and    received    back   at    the   surface   transducer.      The   time   of 
travel   of   the   sound    pulse    is  divided   by   two,    and   this   value 
is    multiplied    by    the    assumed   harmonic    mean    speed   of    sound   in 
sea    water,    thus,    giving   the    distance   according    to   the 
expression   D   =   V    t   .       The   echo   sounder    makes   this    transfor- 
mation  electronically   or   mechanically    within  the   device 
itself   and    displays    the   distance  in   the    water    between   trans- 
ducer   and   transponders.       This    distance    is   not   equal    to    the 
true    distance,    since    the    assumed   harmonic   mean    sound   speed 
generally    does   not  egual    the    true   mean    speed   for   a    partic- 
ular   depth.       To   determine    the    true    harmonic    mean    sound    speed 
one    must    know    the    sound   speed    throughout    the    water    column 
and    apply   it   to   solve   for    the    harmonic    mean.      For    a    distance 
measurement   to   be    accurate,    precise    measurements   of    time   and 
the    harmonic    mean    speed   of    sound   in    the    water    column    are 
mandatory. 

B.  THE  VELOCIT*  OF  SOUND  IN  THE  OCEAN 

The  velocity  of  sound  in  sea  water  is  a  function  of 
temperature,  pressure,  and  salinity:  V  =  V(T,S,P).   Because 
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tne    sound   velocity  is  not    constant    with    position,    sound   rays 
are    refracted  according    to    Snell,s    law. 

We   have   followed    the    standard  practice   of   the   U. S- 
National  ocean   Service    (NOS)    to    use   Wilson's    (1960)     equation 
for    the   speed    of    sound  as   a    function   of    salinity,    tempera- 
ture   and   pressure    (CJmbach,    1976): 


V-1449 . 14+VT+Vp+Vs+VSTp 


(3-1) 


wner  e 


Vm«4 . 5721T-4 . 4532X10~2T2-2 . 6045X10~4T3 

T  -6  4       <3-2> 

+7.9851X10   T^ 


Vp-1 . 60272X10~IP+1 . 0268X10"5P2  {3  _ 

+3.5216X10"9P3-3.3603X10"12P4 


V  =1.39799(S-35)+1.69202X10~3(S-35)2  {3'4) 

o 

VQfTO-(S-35)(-l. 1244X10~2T+7. 7711X10"7T2 

+7. 7016X10~5P-1 . 2943X10?P2+3 . 1580X10"8PT 
+1 . 5790X10~9PT2 )+P( -1 . 8607X10~4T 
+7.4812X10~6T2+4.5283X10~8T3)  (35) 

+P2C-2.5294X10~7T+1. 8563X10"  T2 ) 

+P3(-1.9646X10~10T). 
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In  these  equations  the  absolute  pressure,  P,  is  expressed  in 
kg/cm2,  temperature,  T,  is  in  °C ,  salinity,  S,  is  in  g/kg, 
and  sound  velocity,  V,  is  in  m/s-   For  these  relations  for  V 
the  standard  deviation  from  the  mean  is  0-30  m/s  for  all 
data  obtained  in  the  ranges  -4  <  T  <  30  °C,  1  <  P  <  1000 
kg/cm2,  and  0  <  S  <  37  g/kg  (Wilson,  1960,  p. 1357). 

C.   HARMONIC  MEAH  SOOHD  VELOCITY 

The  slant  range  from  a  ship  to  a  transponder  is  deter- 
mined by  measuring  the  transit  time,  t,  of  an  acoustic  pulse 
and  converting  it  to  distance  using  an  appropriate  value  for 
sound  speed-   If  the  mean  sound  velocity,  V,  in  the  water 
column  is  known,  tne  distance,  L,  from  transducer  to  tran- 
sponder can  be  computed  by: 

L  =  (  t/2  )  V  (3.6) 

For    our    small    project  area    the    distances    between  transpon- 
ders   are   less   tnan  2    nmi,    and   it  is  assumed    that   there   is   no 
horizontal   sound    speed   variation,    and   only   corrections   for 
vertical  sound   speed    variation    are    considered-       For    a 
one-way    travel  time    ,    T,    through  a    water   column    of    depth,    Z, 
the    mean  sound    velocity  is: 

V    =    Z    /    T  (3.7) 

The  meaning  of  "mean  sound  velocity"  is  neither  "mean 
velocity  from  the  surface  to  the  stated  depth"  nor  the  "mean 
value  for  the  velocity  of  sound  through  the  vertical  water 
column".   According  to  Equation  (3-7)  the  mean  velocity  is 
obtained  from  tne  true  depth  divided  by  the  one  way  travel 
time. 
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Figure  3-1    A  Layered  Ocean  there  ?:  =  V^  (z)  Is  the  Sound 
Speed  in  the  i  th    Layer  and  z  Is  Depth- 


Consider  now  a  layered  ocean  having  a  different  sound 
speed  in  each  layer-   Figure  3.1  shows  a  series  of  finite 
layers  of  thicKness  az. ,  each  with  an  associated  sound 
speed,  V-  .   The  time  interval,  A  T.  ,  for  the  sound  wave  to 
pass  vertically  through  the  i  ijj    layer  is  ; 


AT-  =   AZ.  /  V- 


(3-8) 


Summing  tne  time  intervals  for  all  n  layers  in  a  vertical 
water  column  from  surface  transducer  to  bottom  transponder  , 


T  =  IAT.  = 


(3.9) 


Foe    a   continuous   function,    V   =    V(z),    this    becomes 
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I 


\    dz/V(z)  (3.10) 


Equation    (3-7)    may  now  be   written   : 

7 
V    =  - (3.  11) 

'zdz 


I 


,    V(2) 


This    is   the   integral    form   of   the   harmonic   mean-       Using 
Simpson's  Rule,    Equation    (3.10)    can   be    rewritten   as    : 


T    =    (VV„  +4/V,  +2/V2  +4/V3  +     •••    +1/V2)     d/3  (3.12) 

where        d  =   D/n       (magnitude    or    depth    increments) 
D   =   bottom  depth    at    transponder 
n   =   even    number  of    depth   increments   between 

transducer  and    transponder 
V;     =    measured    sound    speed   at   depth    j. 
Since   the    harmonic   mean    sound   velocity    is   obtained   by 
Equation    (3.11),    the    slant    range   from   transducer    to    acoustic 
transponder   can    be   obtained    using   Equations    (3.10)    and 
(3.11): 

L    =    T    V 

CTD  casts  to  140  0  m  were  made  on  18  and  20  May  1985  in 
tne  project  area  just  prior  to  and  during  our  main  experi- 
ment.  Sound  velocities  and  harmonic  means  based  on  the 
sound  velocity  profiles  derived  from  the  CTD  casts  are 
plotted  in  Figures  3.2  and  3. 3  .   Although  the  profiles 
shown  in  these  figures  are  for  different  days  and  different 
places  within  the  project  area,  they  shoh  very  little  varia- 
tion in  harmonic  mean  sound  velocity,  and  thus  support  our 
previous  assumption  that  the  sea  water  is  horizontally  homo- 
geneous to  an  acceptable  degree  within  the  study  area. 
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Harmonic   Mean   and   Sound    Velocity    Profile   2. 

26 


HARMONIC  MEAN  COMPARISON 


O     H82.97 


K83.01 


146  3.05 

r  i  . 


K83.09 


B 


o 

3- 


z2 


18   MAY   1985 


20   MAT    1985 


i         i         i 


HS2.97 


Kro.Ol  H63.05 

HARMONIC   MEAN   VELOCITY    CM/5) 


— i 

1453.09 


1 


Figure   3- 4        Harnonic   Mean    Comparison. 
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Actually,  the  difference  of  these  two  harmonic  mean  veloci- 
ties is  approximately  0.07  m/s  near  the  bottom  as  is  shown 
in  Figure  3.4.   The  mean  of  these  two  harmonic  mean  sound 
velocity  profiles  is  used  in  the  range  calculations  of  the 
following  chapter. 

D.   ACOUSTIC  WAVE  HEFBACTION 

To  simplify  our  acoustic  analysis  it  is  desirable  to 
assume  that  in  the  study  area  sound  follows  straight  paths 
from  the  surface  to  the  ocean  bottom.   To  test  this  approxi- 
mation the  differences  in  path  lengths  from  the  surface  to 
the  bottom  were  calculated  for  varies  depression  angles  at 
the  surface. 

Frofiles  such  as  those  of  Figures  3.2  and  3m3   are 
usually  simplified  for  analysis  by  separation  into  an  appro- 
priate number  of  segments  each  having  an  approximately 
constant  gradient  (Figure  3.5).   If  the  velocity  of  sound 
changes  linearly  with  depth,  sound  rays  can  be  shown  to  nave 
a  constant  radius  of  curvature  (UricJc,  1983,  p.  124). 
Horizontal  distance,  d,  depth,  z,  radius  of  curvature,  R, 
segments  of  arc,  s,  and  depression  angle,  A,  are  shown  sche- 
matically in  Figure  3.6.   The  path  of  a  ray  when  the  speed 
of  sound  varies  witn  depth  can  be  calculated  by    application 
of  Sneli's  law: 

cosA,  /  C,  =  cosA2  /  C2  =  •••  =  cosAV  C-  (3.13) 

wnere   C^=  sound  speed 

Ar  =  depression  angle 

For  an  initial  depression  angle  of  30°  the  profile  of 
Figure  3.5  gives  the  results  shown  in  Table  I. 

The  raaius  of  curvature  can  be  obtained  in  the  following 
manner  (Kinsler,  1982,  pp399-402).   Let 
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Figure  3.5    Simplified  Sound  Velocity  Profile- 
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Figure  3.6    Schematic  Diagram  of  Sound  Ray  Trace. 
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TABLE   I 
Depression  Angle   at  Different  Layers 


I 

Sound           | 

Depth 

1      Depression 

Points 

I 

Velocity         I 

(m) 

|    angle    (deg.) 

1- 

I 

1493.0           | 

0- 

I               30.  0 

2. 

I 

1483.5            j 

176.4 

|               30.63 

3. 

I 

1484.04          | 

266.  6 

I               30.59 

4. 

I 

1481.2            | 

526.4 

|               30.78 

5- 

I 

1482.46         | 

1097.2 

i              30.69 

6. 

I 

1484.34         j 

1400.0 

I              30.57 

J 


9":    = 


(    C.    .   -    C     )    /    {    z.       -    2  .    ) 


(3.14) 


and 


E;    =   -     (C.    /    g.  )     cosA- 


(3.15) 


where 

fi-    =  Radius   of    curvature   in    the   i*ij    depth    increment 

C;     =   Sound   speed   in    the   itjj    depth    region 

A^     =   Depression    angle 

g.     =   Speed   gradient 
The    change    in   range   Ad.    and    depth  4z-    are    (Kinsler,     1982) 


AdL    =   -hl    (    sinA-    -   sidA.+)  ) 


(3.16) 


AZ-      =     -E;       (     COS  A;+      -     COSA-      ) 
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(3.  17) 


so    the   chord   can    be   obtained      by   two    components    of    change  in 
range   and   depth;       that  is, 

D.     =    (  Ad.2    +AZ.2    p/2  (3.18) 

The    angle,    0-  ,    is 

9-     =    2sin~i     (   D-  /2R .    )  (3.19) 

fa  J.  C 

so    the   arc,    s.,    is: 


L 


s-    =    R-  9-  (3.20) 

where      0;     is   the    angle   in   radians. 

The    range   and    depth    increments    give   the    straight   line 

L    =   [  (2AdL    )2+(2^z-     )  2]1/2  (3.21) 

The  ray  paths  can  be  obtained  by  summing  arcs: 

S  =  ls-L  (3.22) 

The      difference  due      to      refraction      between   a      sound      ray's 
actual   path    and    a   straight    line   is    thus: 

E    =   S    -    L  (3.23) 

where  E    =   the    error  between   two   lines, 

S    =    the    distance   along    the   ray    path,    and 
L    =    the    straight-line    distance    between    the 

ship-mounted    transducer   and    bottom    transponder- 
Ray-trace   computations    were  made   for    various  initial 
depression    angles    (at    the    surface)     down    to   a    lepth    of    1400    m 
(Tanles    II    and   III)Since    E<1    cm    for    a    30°    surface 
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TABLE  II 
The  Result  of  Ray  Trace  Computation 

For  Initial  Depression  Angle  =  30  Degrees 

Depth   Gradient   Radius    Horiz.    Depth     Arc 
Increment   g        R     Dist- (d  )    z        s 

Cm] On] Do] [ml 

1.  0.0538  32011.36    301-72  176.4  349.5023 

2.  0.0060  287966.49    152.47  90.2  177.1543 

3.  0.0109  157706..64    437.35  259.3  509.1273 

4.  0.0022  780984-82    960.03  570.8  1116.8993 
5-  0.0062  277668-87    511.35  302.8  594.2820 

Total    :  2363.42         1400.0       2746.9658 

The   Distance  "of   Straight    Line      :  2746.9562 

The    Error  :  0.0096    meters 


« 


depression  angle  and  a  1400-m  water  depth,  corresponding  to 
a  horizontal  distance  of  2363  m,  which  covers  the  project 
area,  and  since  the  depression  angles  encountered  during  the 
project  were  30°  or  greater,  1  cm  represents  an  upper  limit 
to  the  error  resulting  from  our  assumption  of  linear  ray 
paths  from  the  surface  to  the  bottom,  and  a  correction  for 
refraction  is  unnecessary. 
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TABLE  III 
The  Error  Coaparation  at  Different  Angles 



Depth   :   1400  meters 


Depression 

Horizantal 

Ray 

Straight 

Error 

Angle 

Distance 

Trace 

Line 

(degree) 

(m) 

(m) 

(m) 

(cm) 

0 

13429-23 

13508.3850 

13502.0064 

637.86 

10 

t>637. 03 

6783.2984 

6783.0795 

21.89 

20 

3647.  10 

3906.6140 

3906.5806 

3.34 

30 

2363.42 

2  74  6.96  58 

2746.9562 

0.  96 

60 

801.28 

1613.0859 

1613.0853 

0.06 

90 

0.0 

1400.0 

1400.0 

0.0 
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17.  RELATIVE  POSITION  DETERMINATION  OF  THE  TRANSPONDER  ARRAY 

A.  INTRODUCTION 

Two  methods  to  determine  relative  positions  are 
discussed  in  this  chapter.   They  are  basic  baseline  crossing 
and  the  modified  baseline  crossing  method;  the  former  is 
described  by  Hart  (1967)  and  Haehnle  (1967),  the  latter  is 
described  by  Fubara  and  Mourad  (1972).   The  main  difference 
between  the  methods  is  that  Hart  and  Haehnle  require  accu- 
rate transponder  depths,  but  Fubara  and  Mourad-  do  not. 

The  cloverleaf  method  used  to  determine  transponder 
depths  is  described  by  Haehnle  (1967)  is  introduced  in  the 
next  section. 

B.  CLOVERLEAF  METHOD  FOR  TRANSPONDER  DEPTH  DETERMINATION 

An   initial  approximation   to   the   depths   of   bottom-mounted 
acoustic   transponders  can   be   obtained   from   a   knowledge   of 
the    water   depth   at  the   point   where   the    mooring    is  released 
from    the   ship.      However,    the   unit    may    move   laterally   while 
descending,    or   bottom   topography   may   have   considerable 
relief    not    evident  on    the   echo   sounder    record-       Hence,    a 
more    accurate    method    of   depth    determination   is    required- 
Here    the   depth   is    defined    as    the   distance    between    the    hull- 
mounted    transducer   and  the    bottom-anchored    transponder    when 
the    vessel   is    positioned   exactly   over    the    trans ponder- 

The    ship   need   not   be    directly   above   the   transponder   but 
should   be   in   close   proximity    to  it.       Although   it    is    very 
difficult   to  cross  over   the    transponder   exactly    with    the 
ship,    the   depth   may    be  approximated  as    the    minimum   slant 
range   when    the   transponder    is   crossed   over    many    times   using 
a   "cloverleaf"    method    (Figure    4.1).       Thus,    if    a    transponder 
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b.      Maneuver  for  obtaining   four  CPA's 


Figure    4-1         Deteraination    of    Transponder    Depth. 
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D=50  m 


CPA 


Figure   4.2         Exauple   of   Depth   Difference. 

is    at    a    depth    of    153  0    m   and    the   ship   is   displaced   50    meters 
horizontally   from   it,    tne    depth    determined   by   the   cloverieaf 
method    would   be   in   error    by    only   0.33    meters    (Figure    4.2). 

The    procedure   used   for    running    a   cloverieaf    maneuver    is 
as    follows    (Haehnle,     1967)  : 

1.  Each    time    a   transponder    is    dropped    from    the   ship   the 
ship's   position  is   noted. 

2.  The    ship    proceeds    toward   the    transponder    drop    point 
and    passes    near    it;       A    closest    point    of    approach 
(CPA)     is   then    determined   and   recorded. 

3-       The   ship   then    turns,    crossing    the    first    track    perpen- 
dicularly   and    passing    near   the    first   CPA.      Once 
again,    the   CPA   of    this    run   is    recorded. 

4.       Tne    above    maneuver    is    repeated    until   at    least    four 
such    CPA's    have  been    obtained,    each    time    attempting 
to   pass   through   the    previous  CPA    perpendicular   to    the 
track. 
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Figure   4.3        Saiple  of   Depth    Determination. 

5.      The   depth    of    the   transponder   is   determined    by    plot- 
ting   the   slant    ranges    versus    time    for   each    of    the 
four   CPA   runs.      The    minimum    point    of   t^is  curve   is 
considered    as    the    minimum   slant    range    for    that 
particular    transponder     (Figure    4.3). 
The    R/V    ACA:JIA   ran  similar    cloverieaf    patterns    on    21    Say 
1985.      Sne    crossed   near   each   transponder   many    times    as 
closely    as   possible.       The    minimum    distances    and    corrected 
depths    are    listed    in    Table    IV    . 


38 


TABLE   17 
Transponder  Depth    by   Cloverleaf   Maneuver 

(From    Oceano   [1985,    p-12]) 

Corrected 
Transponder  Depth 

(m) 

1  1366.9 

2  1331.3 

3  1381.5 

4  1406.7 


j 


i 


C.        BASELIHE    CfiOSSING    METHODS 

Baseline    length    can  be    determined    utilizing    the    baseline 
crossing   method   described   by   Haehnle    (1967)    and    Hart  (1967)  r 
which   involves   repeatedly   transiting   the   baseline    while 
simultaneously   measuring    ranges   to    two    transponders.       The 
minimum   sum   of    the   two  horizontal   ranges   as   determined   from 
measured   slant   ranges   and   known    transponder   depths   is   taken 
to    be    the   baseline   length. 

The    basic    baseline  crossing    method    requires    a   knowledge 
of: 

1.  The   ship*s  course    and   speed, 

2.  The   accoustic   slant   range    (or    the   two-way    travel 
times   from    the    ship    to   two   adjoining    ocean- 
bottom   transponders,    and 

3.  The    harmonic    mean    sound    velocity. 
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Figure  4.4    Saup 


le  of  the  Baseline  Crossing. 
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The   ship,    travelling   at    a   constant   speed  and  heading, 
crosses    the   lines   joining    the    transponders  as   shown   in 
Figure   4.4.      Tha    acoustic   slant   ranges    are   recorded    as    func- 
tions  of   time    during    the   crossings.      The    sum   of    two   hori- 
zontal  ranges    (the  distance    T1-T2    in   Figure   4.4)    is    minimum 
when    the  ship   is    in    the   vertical   plane   containing    the   two 
transponders.      Similarly,    for   a  four- transponder  array,    the 
distances   T1-T2,    T2-T3,    T3-T4,    T1-T4,    T1-T3    and    T2-T4    are 
determined.      The    minimum  sums   can    be   obtained   by   least- 
sguares   fitting   the    horizontal    ranges    against    time   and 
obtaining  the   vertex    of    the    parabola   as    shown   in   Figure    4.5 
During   the    experiment   four    to   six   crossings    were   made    for 
each    baseline    between    transponders.       The    results   are 
discussed   below  in  Sections    1    and   2. 

Although   the    baseline   crossing    method    is  relatively 
simple   with   respect    to   data   reduction,    it   is   operationally 
time    consuming.      During   the    experiment   it   took    one    hour    and 
thirty   minutes  for   each   baseline   crossing.      Thus   the   eight 
crossings    (4   for   north-south   crossing,    4    for   east-west 
crossing)     took   a    total  of    12    hours. 

1  -      i&sic    Baseline   Crossing   Model 

Figure   4.4   shows   a    ship   track    crossing    a   baseline. 
The    ship's   transducer   was   set    to  interrogate   the   bottom 
acoustic   transponders  every    10    seconds   for   this    experiment. 
Ranges   were   obtained    from   each    transponder-      It    is   easily 
seen     (Figure   4.4)     that,    as    the   ship   approaches    the    oaseline 
I1T2,    the   total   horzontal   range  at    the    first   ship   point    (0) 
is    greater    than   at   the  second    ship    point    (1),    and   so   on, 
until  a    minimum  sum    is  reached   at   the   instant   of   baseline 
crossing;    from   there    on   the    total   increases   for    each    point 
on    tne   track.      Since    the    ship    maintains    a   constant    course 
and    speed,    the    mathematical   function    which    describes    tnis 
relationship   is   a    parabola    (Figure    4.5)    of    the    form     (Hart, 
1967)  i 

4  1 
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Figure   4.5 


Sample    of    Range-Tine    Parabola    Showing 
baseline    Crossing   Time. 


H=AT2    +    BT    +    C 


(4-D 


wnere 

A,  B,  Z    are  constants  to  be  determined, 

T  is  time  [  s  ],  and 

E  is  total  horizontal  range  between  transponders 

To  determine  the  unknown  variables  A,  B  and  c  in 
Equation  (4.1)  a  least  squares  adjustment  is   performed.   At 
least  four  sets  of  data  R  =  R(T)  are  needed  to  solve 
Equation  (4.1).   Since  observed  parameters  contain  random 
error.  Equation  (4.1)  can  be  rewritten  as: 
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A   12    +    B    T    ♦    C    =   H    +    V 


(4.2) 


where 


V   is   the   residual   of    the    observation. 


In    matrix   notation   this   corresponds   to, 


*z 


+ 


-T/,-T2^/ 


which   may   be   written    as 


> 

/       \ 

f*,] 

! 

A 

Rz 

8 

= 

) 

c 

■ 

^  Rr»< 

(4.3) 


V         +    T         X         =   ft 
TIM        /ix3    3*/        nxi 


(4.4) 


where        V-     are   the  residuals   of   the   ooservations, 
T.  .     are  the    numerical    coefficients   of   the 

unknown   variable, 
X:    are  the   unknown    variables (A,    B,    C)  ,    and 

d 

R.    are  constant    terms    of   total   horizontal   range. 

The   standard    deviation    of   the   mean   of   sound    velocity 
is    ±0.3    m/s    (Wilson,     1960).       So,    the   standard   deviation    of 
the    distance    is    the    following,    where    time    is  assumed    error- 
less   : 


<k.  -   tfc     I 


(4.5) 


Each  observation  is  considered  to  be  independent  and 
uncorrelated.   Since  a  measurement  of  high  precision  has  a 
small  variance,  and  one  of  low  precision  has  a  large  vari- 
ance, the  higher  the  weight  the  higher  is  the  precision  and 
vice  versa.   Accordingly,  the  weight  W  of  a  single  obsera- 
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tion    is    defined  as  a    quantity   that   is   inversely    proportional 
to   the   variances    of    the   observations    (d2)     (Mikhail,    1981, 
p.  6  6)  . 


W    =  d2    /  6  * 
o 


(4.6) 


where    d2  =  variance  of  unit  weight 

o  3 

For    the      adjustment    of      Equation    (4.4),         the   least      squares 
criterion  requires      minimization   of   the    weighted      function    G 
(Mikhail,    1981,    pp.69-73): 


=    W 


V,2    +W    V2     ♦    •••    +W    V2    =     1  F.  V.2 

'    '         z   z  ft  n.        is/    t    l 


(4.7) 


which   in   matrix  form    is   written   as 


G  =  V   »  V 


(4.8) 


where 


W  - 


X  0  •  •  0 


lA/, 


^  ; 


n*rv 


Now  Equation  (4.4)  can  be  rearranged  as 


V  =  R  -  T  X 


(4.9) 


Substituting    for    V   into   Equation    (4.8),    get 


t  t      t    ± 

G    =     (R-TX)     W(R-TX)     =     (R    -X    T    )  H  (R-TX) 


(4.  10) 


so, 


f  t     t  i-  -k    + 

3    =    R    fc/R-X    T    Wfi-B    fcTX+X    T    WTX 


(4.  11) 
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Since  G    is   a   scalar,       the   right   hand   side   of    Equation    (4.10) 
is   also   a  scalar.      Furthermore,    the   transpose   of   a   scalar   is 
equal   to   the   scalar: 

XtTtWR    =     (XT    WE)     =    R    WTX  (4.12) 


and 


W      =    W  (4.  13) 


Thus,    Eguation    (4.10)    can   be    written   as: 

G    =    E*WR    -    2RtWTX    ♦    X     (I    II)  I  (4.14) 


In  Equation  (4.14)  all  matrices  are  constants,  except  X,  the 
matrix  of  unknowns.  For  G  to  ce  a  minimum,  its  partial 
derivative    with   respect   to    X   must   be   equal   to   zero: 

£&  „  t  „    t   «t 

=    -2B    WT    + 


ax 

and , 


2X    (T    WT)     =0  (4.  15) 


t  t 

(",„„  W^J         )     X         =    T         W  R  (4.  16) 

3xn    I**   n*3  jxi  l*n    run    n-X/ 

giving   the   solution    for   X: 

t  t 

X    =    (T       W    T    )  -»    T       W    R  (4.  17) 


Now,  the  minimum  distance  occurs  when  dR/dT  =  0 
which  is  given  below: 

T  =   -  B  /  2  A  (4.  18) 
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Thus,    to  get   the    minimum   distance   of    the    sum,    one   may 
substitute    the   time   from    Equation    (4.18)     into   Equation 
(4.1)  : 

H     (min.)    =   B2/4A    -    B2/2A    +    C  (4.19) 

where        A,    B,    C   are    given   by   Eguation    (4.17). 

The   variances   of    A,    B    and  C    (Mikhail,    1981,    p258) 
are: 

di2    =   6*    (T   HI)-»  (4-20) 

The    standara   deviation  of    T   in    Equation     (4.18)     is: 

dT     =    [  (dA    B/2A*)2    +    (tfB/2A)2]*/>  (4.21) 

The    standard    deviation  of    the    minimum    distance    R    in   Equation 
(4.  19)    is: 


6    =[  (T2Cf   )2+  (tu    )2-k52+  (2AT*B)2($2  -ji/2  (4.22) 

R  A  q  C 


2-      Data   Processing 

During   the  experiment    at   sea,    both    north-soutn    and 
east-west   baseline  crossings   were   run    (Figures    4.6    and   4.7) 
using    data   collected    by  Oceano   Instruments    (1985a). 
For    each   baseline   crossing,     11    consecutive    interrogation 
points   were   fit    to   a    parabola    by   least    squares.      The   reason 
foe    using    11    points    is  that    this   number    is   sufficient   to 
define    the    parabola    near    its    minimum.       The    minimum    value   of  the 
parabola   can    be   considered    to    define    tne    baseline    length. 
The    parabolas    shown    in   the    Figures    4.8,    4.9,    4.10,    4.11, 
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Figure   4.6         north-South   baseline    Crossing. 

4.12    and   4.13     (at    the    end    of    this   section)     are    for    single 
crossings    between    X1-T2,    12-12,    T3-T4,    T1-T4,    T1-T3    and 
T2-T4.      Tne    results    are    tabulated   in    detail    in    Table  V,  which 
gives  tne   baseline   lengths  and    their    standard    deviations   ior 
each    crossing.       To   get   the    mean   of    these    lengths,    one   can 
use    the    following    equation: 


ft       =    (E1+E2*    •••    ffin)     /   n    =    f  (fi) 


(4.23) 


Since    the   standard    deviation    of    the    mean  length    ft   in 


m 


Equation    (4.23)    is: 

6m  =[  (c5i  a/c/9Ei)2+  (62  af/aR2)2+  •••    ]l/2 


(4.24) 


we    nave 


47 


Figure   4.7        East-Iest   Baseline    Crossing. 


0^    =[  (6l/n)  2  +  (d2/n)  2*»««+  (tfn/n)  2  ]l/2 


(4.25) 


Equations  4.23  and  4.25  are  used  with  the  data  of 
Table  V  to  compute  the  mean  of  baseline  lengths  and  their 
standard   deviations         (Table    VI)» 

lyman    BurJce    of   Oceano    Instruments    and   L.    Spielvogel 
of    Seaco,    Inc.,    have    also  computed    the    baseline    lengtns. 
They    use   a    different    method    and   technique    using    most    of    tne 
data    collected   during   the    experiment.      Their   results   and   the 
differences   between    them    are   tabulated    in   Table    VII-    Due   to 
our    more   limited    data    set   and    uncertainty   in   the   depth    for 
the    baseline   crossing    method,    these   results    differ    from 
tnose    obtained    with    the    crossing    method    by   about    2    ra. 
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TABLE    7 

Baseline  Lengths   Determined   by  Baseline  Crossings    (m) 

Crossing  J             T1    -    T2  |             T2    -    T3  |  13   -    T4 

1-  |     1793.32±2.43  \  1993.09±0.61  |  1972-39±2.50 

2.  |  1792.62±2.17  |  1991.92±2.  73  |  1969.50±2.02 

3.  |  1793.67±2.27  |  1992.75±2.38  |  197  Q-29±2.  1  9 

4.  |  1791.77±2.22  |  1994.33±2.09  |  1969. 97±2. 23 

Crossing  j            II    -.14  |           11-13  J  12  -    T4 

1.  |  1860.8U3.25  j  2789.62±2.20  |  2588.66±2.03 

2.  j  1858.43±1.81  |  2789.77±1.62  |  2589.70±1.94 

3.  i  1858.63±2.23  I  2789. 45±2.  11  |  2587.88±1.14 

4.  |  1858.48±2.08  |  2789.73±1.63  |  2588.17±2.21 

5.  |  1858.65±2.04  |  2790.22±2.26  |  2589.22±2.2b 

6.  |  |  |  2588-82±2.20 


i 


3 .   Modified  Baseline  Crossing  Method 

Mourad  and  Fubara  (1972)  developed  this  method. 
Their  method  determines  the  relative  positions  of  ocean- 
bottom  transponders  and  the  depth  of  each  trarsponder- 
There  is  no  necessity  to  accurately  measure  each  depth 
directly;  this  is  the  main  difference  from  the  previous 
work.   The  physical  principle  of  this  method  is  shown  in 
Figure  4.14   S  is  the  baseline  length  and  L1  and  L2  are  a 
pair  of  coplanar  ranges  from  the  ship  to  transponders  1  and 
2  whose  depths  are  Z1  and  Z2,  respectively.   The  mathmatical 
expression  for  this  configuration  is; 
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r   ■■  ■ 

■  ■  ■    -  -               -  -- -■■  ■       i 

I 

TABLE  71                           1 
dean  of  The  Baseline  Lengths             I 

The 

1    Mean  Value    1   Standard  Deviation 

1      (  m  )       1         (  ffl  ) 

"  "  ""  "■■'    " j 

T1  - 

T2 

|     1792.84      j        1-14 

T2  - 

T3 

|     1992.42      |        0.93 

T3  - 

T4 

|     1970.54      |        1.12 

T1  - 

T4 

|     1859-00      i        1.04 

T1  - 

T3 

|     2789.76      |        0.89 

T2  - 

T4 

|     2588.74      J        0.82 

L 

j 

S    =     (    L12    -    Z12    )i/2    +     (    L22    -    Z22    )W 


(4.26) 


To  determine  the  three  unknowns  S,  Z1  and  Z2  by 
least  squares  method,  at  least  4  pairs  of  Li  and  L2  are 
needed-      Now,    Eguation    (4.26)    can  be    rewritten    as   follows: 


F    =    S- (    L12    -    Z12    )  1/2    -     (    i22    -    Z22    )i/2 


(4.27) 


Assume 


A1     =     (    L12    -    Zl2    )   1/2 


(4.28) 


A2    =     (    L22    -    Z22    )  1/2 


(4.29) 


Equation    (4.27)     is   linearized    and    partially    dif terentiated: 
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TABLE    YII 

Baseline  Length      Determination   by   Oceano 
Instruments   and   L-    Spielvogel 


Lines 

Distance         (m) 

Difference 

l             Oceano 

I        Spielvog 

el 

I            (m) 

T1    -   T2 

J            1792.68 

I            1791.7 

|            0.98 

T2    -   T3 

|             199  7.46 

|            1995.9 

|            1.56 

T3    -   T4 

|             1973.92 

I            1972.6 

I    -        1.32 

T4    -   T1 

|             1861.86 

|             1860.5 

1.36 

T1    -   T3 

j            2792-69 

|            2790.8 

|            1.89 

T2    -   14 

I            2592.18 

I            2590.3 

1.  88 

2L 


=    -A1-1L1 


(4.30) 


=    -A2-1L2 


(4.31) 


9P 
3S 


=    1 


(4.32) 


221 


=    A1-1Z1 


(4.33) 


9Z2 


=    A2-1Z2 


(4.34) 
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Figure   4.8        The   Best   Fit    of    Parabola    for   T1-T2  Crossing. 
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T2-T3  BASELINE  CROSSING 


20 


30 


"i — ' — r 


T 


i — ' — r 


T 


1 


40  50  60  70  90 

TIME  IN  SECONDS 


90  TOO         110         T20 


Figure    4.9        The   Best  Fit   of    Parabola    for   T2-T3   Crossing. 
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Figure    4.10        The   Best  Fit    of    Parabola    for   T3-T4   Crossing. 
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Figure    4.11         The   Best  Fit    of    Parabola    for    T1-T4   Crossing. 
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Figure  4.  12    The  Best  Fit  of  Parabola  for  T1-T3  Crossing. 
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Figure    4.13         The   Best   Fit    of    Parabola    for    T2-TU    Crossing. 
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L. 


Figure   4-14        Baseline   Crossing   Configuration. 

Tne    resultant    ooservation    equation    for   a    least- squares    solu- 
tion   is     (ilourad    ,     1972,    p. 21): 


AA    +    BV+W    =    0 


(4.35) 


which    can    be    written    as. 


f  /    A\''z\.  Ai'/Z 


a5 
aZZ 

•  jl  - 


(      -I 


-AlLl,f~A2tf,0 

O    .     O   -AUI2 


(w?    f  ^ 


v\. 


•  / 


K  vv» 


ru 


0 


sO 


(4.36) 


Tnus,    tne   xeas t- squa res   solution    (Mourad,    1972,    p. 21)    gives: 
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A  ^[P^+A*  (BP-iB*) -*&]-*&*  (BP-»Bt)-*H  (4.37) 


V    =-P-iBt  (BP-iB*)-!  (AA4-W)  (4.38) 

The    variance   of   unit    weight    d0  isi 

6Q  =    {[(BP-*Bt)-i    (AA+W)  ]tW/df}  1/2  (4.39) 

The    variance-covariance   matrix    for    the   adjusted    parameters    Xa 
is: 

(£  £P^  +A*  (BP-iBt)-»A]"i  (4.40) 

The    adjusted    parameters   X^  is: 

Xa    =    X0    +  A  (4.41) 

where 

A  is  a  coefficient  matrix  for  the  unknowns  correction, 
B  is  a  coefficient  matrix  for  the  observation, 
v  is  a  vector  of  residuals,  representing  the 
corrections  to  observed  ranges,  L1  and  L2 , 
A  is  the  correction  to  assumed  X0 , 
W  is  a  constant  matrix  in  Eguation  (4.27),  when 
using  assumed  value  X0  (S0  ,L10,L2o), 
P(  is  the  weignt  matrix  of  L0, 

Px  is  the  weighting  function  associated  witn  XQ , 
60    is  tne  variance  of  unit  weight 
Xa.  is  true  value  of  S,    Z1  and  Z2 ,  and 
X0  is  an  approximate  value  of  Xa. 

The  weignt  matrices  are  taken  to  be  unity  as  the 
measurements  were  of  egual  precision.   Tanle  VIII  shows  the 
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resulting  horizontal  distances  between  transponders  and  the 
depths  of  each  transponder.   These  results  are  obtained  by 
measuring  pairs  of  slant  ranges  at  the  instant  the  ship 
crosses  the  baseline.   Since  the  ship  crossed  the  baseline 
four  times  for  the  baseline  crossing  method,  there  are  only 
four  pairs  of  slant  ranges  which  can  be  used;  these  repre- 
sent the  smallest  data  sets  needed  for  a  least  squares  solu- 
tion. 

The  baseline  lengths  computed  by  Oceano  Instruments 
and  Dr.  Spielvogel  are  given  in  Table  IX   The  difference 
between  Oceano's  and  Spielvogel's  computed  values  are  witnin 
the  precision  of  the  system.   My  results  in  Table  71. II 
differ  from  Table  IX, since  Oceano  and  Spielvogel  used  much 
larger  data  sets  (more  than  250  data  points) ,  whereas  I  used 
tne  much  smaller  sets  (11  points  for  each  of  four  crossings) 
required  by  the  two  baseline  crossing  methods,  namely  eleven 
data  points  for  each  of  four  crossings.   Since  Oceano  and 
Spielvogel  used  a  resection  method  with  four  transponders, 
it  is  not  necessary  for  tne  ship  to  travel  the  naseline 
repeatedly:   The  only  requirement  is  that  the  ship  De  in  the 
vicinity  of  the  center  of  the  experimental  area,  where  data 
can  be  collected  continuously.   On  the  contrary,  the 
crossing  methods  require  a  large  amount  of  ship  time  to 
collect  a  small  amount  of  data.   Thus,  it  is  not  surprising 
that  there  should  be  good  agreement  between  Oceano  and 
Spielvogel  and  rather  poorer  agreement  between  my  results 
and  theirs. 
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TABLE  Fill 
Results  of  Modified  Baseline  Crossing  Method 


Trans- 

Horizontal          | 

Transponder 

Mean 

ponder 

Distance    Between    J 

Number- 

Depths    (  m) 

Pairs 

Transponders    (m) | 

1-2 

1782.91±9.51              | 

1 

1361.  52±4-53 

2-3 

1993.02±6.30              | 

2 

1333- 94±2.54 

3-4 

1964.34±8.32              | 

3 

1375- 76±5.21 

1-4 

1836.65±9-21              j 

4 

1407. 08±3-47 

1-3 

2815.47±15-60            \ 

2-4 

2574.13±7.55              I 

I 
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TABLE  IX 

Comparison  of  Results  for  Different  Methods 

Basic  Modified 

Lines   Oceano    Spielvogel    Baseline  Baseline 

Crossing  Crossing 

1792.84±1.14  1782.91±9.51 

1992.42±0-93  1993-02±6.30 

1970.54±1.12  1964.34±8-32 

1859.00±1.Q4  1836.65±9.21 

2789.76±0.89  2815.47±15.6 

2588.74±0.82  2574.13±7.55 


T1-T2 

1792-68 

1791.7 

T2-T3 

1997.46 

1995.9 

X3-T4 

1973.92 

1972.6 

T  4-T  1 

1861.86 

1860.5 

T1-T3 

2792-69 

2790.8 

T2-T4 

2592.18 

2590.3 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 

Two  methods  were  used  to  determine  the  horizontal 
distances  between  the  elements  of  an  array    of  four  acoustic 
transponders  lying  on  the  ocean  bottom  at  an  average  depth 
of  1370  meters,  namely  basic  and  modified  baseline  crossing 
methods.   The  principal  differance  between  the  two  tech- 
niques is  that  the  former  requires  precise  depth  informa- 
tion, whereas  the  latter  does  not.   A  "cloverleaf"  maneuver 
was  described  which  may  be  used  to  determine  the  depths  of 
the  transponders.   Because  a  larger  number  of  acoustic  range 
observations  was  used  for  the  basic  method,  the  standard 
deviations  for  the  horizontal  distances  between  transponders 
using  that  techniaue  were  about  one-tenth  those  found  using 
the  modified  method.   Our  results  for  the  basic  method  are 
consistent  with  independent  calculations  made  by  Spielvogel, 
who  used  a  much  larger  data  set.   The  differences  between 
our  results     using  the  modified  method  and  Spielvogel's 
are  within  our  relatively  large  standard  deviations,  with  the 
exception  of  the  distance  measured  between  transponders  1 
and  3.   This  discrepancy  is  probably  due  to  the  ve.ry    oblique 
crossings  of  the  1-3  baseline  which  were  used. 

In  addition  to  the  smallness  of  the  data  sets  used  for 
our  calculations  other  sources  of  error  in  our  results  may 
stem  from  the  variability  of  the  ship's  speed  and  its  pitch 
and  roll,  which  we  did  not  take  into  account. 

If  the  absolute  positions  of  the  bottom  transponders  are 
required,  the  position  of  the  ship  may  be  determined  by  an 
absolute  positioning  system  such  as  GPS  while  measuring 
slant  ranges  from  the  ship  to  each  transponder.   When  at 
least  three  or  more  absolute  positions  of  the  ship's 
acoustic  transducer  are  known,  the  transponder  coordinates 
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may  then  be  adjusted  by  the  method  of  least  squares  to 
obtain  a  best  fit  between  measured  and  calculated  slant 
ranges  from  the  known  transducer  positions.   This  method  is 
introduced  in  the  Appendix. 
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APPENDIX    A 
ABSOLUTE    POSITION    OF    OCEAN    BOTTOM    TBANSPONDERS 

The    absolute    position   of   a   ship   can    be   obtained    from   the 
Global  Positioning  System    (GPS).      Here,    ship   position  means 
the    position   of   the    GPS   antenna-      The   antenna   position   is   a 
function  of   time-      Once  the    absolute    position   of   the   antenna 
is   obtained,    and    simultaneously,    the   ship's   spaed,    heading, 
pitch,    roll  and   acoustic   data,    then   the   absolute   position   of 
antenna    can   be   used    to  calculate   the   absolute   position   of 
the    bottom    transponders-      The   procedures    are: 

1-      To   determine    the  absolute   position   of    the   GPS 
antenna; 

2.  To  convert    the   absolute    position   of   antenna    to   the 
hull-mounted    transducer    using   the    ship   heading,    pitch 
and   roll   data ;    and 

3.  Simultaneously,    to   convert   the   ship   transducer's 
position   to  the  oottom    transponders  by   using    acoustic 
data. 

A   method   is  introduced    here   to   determine   the   absolute 
position  of   sea-floor    transponders,    if    the   absolute   position 
of    surface   transducer   is   given-      Least-squares    adjustments 
must    also  be    supplied   to    get    the  best    fit   of   the   transponaer 
posit  ion. 

The    three-dimensional   geometry    used    is   iilustrated   in 
Figure   A.1-      There   are   4    transponders    (T1,    T2,    T3  and   T4)  . 
Their  estimated  coordinates   are    (X-,Y^#Z*),    the    absolute 

position   of   ship-mounted    transducer   is   S  (0«  « V-  #  W.  )  ,    and   tne 

J       J       J 

four    slant   ranges   are   i». 

la 

The    adjusted    distances    %     between    transducer   and    tran- 
sponders   can   be   computed   given    coordinates   of    any   two 
points: 
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:x?y^z3) 


Figure   A. 1 


The   Belationship   of    Transducer    (   S   ) 
To    Transponders,     (   T«).« 
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2(    =   [  <{J-X1)2+(V-I1)2«-  (H-Z1)  2  ]*/2  (A-  1) 


A. 


J?     =  [  (0-X2)2+(v-Y2)2+  (»-Z2)2]i/2  (A. 2) 


£.    =   [  (U-X3)  2+(V-Z3)2+  (K-Z3)  2]i/2  (A. 3) 


3 


y"x 


J      =   [  (0-X4)  2+(V-Y4)  2+  (W-Z4)  2]i/2  (A. 4) 

which   can   be   written    simply    as    follows: 

\    =   n0j-XL)2+  (V. -Y.  )2+(w.  -z-  )2]i/2  (A.5y 

wnere   i  are  1,2,3,4.   index  corresponding  to 
each  transponder; 
j  are  1,2  *••  N.   transducer  positions 
at  the  surface; 
U# V, W  are  coordinate  of  transducer; 
X,Y,Z  are  coordinate  of  transponder;  and 
i.  are  adjusted  distance  of  slant  range. 

It  was  necessary  to  determine  the  minimum  number  of 
positions  (N)  of  the  hull-mounted  transducer  that  were 
needed  to  fix  the  position  of  four  transponders  in  the 
taree-dimensional  coordinate  system.   There  are  three' 
unknown  variables  to  be  solved  for  each  transducer  position, 
and  there  are  four   transducer  positions  to  solve  for.   To 
solve  for  tnese  twelve  unknown  variables  and  apply  least- 
squares  adjustment,  the  number  of  equations  must  be  greater 
than  cr  equal  to  the  number  of  unknown  variables.   Tuis  is 
given  by; 
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4N    >    12  (A.  6) 

where      4  is   four    equations    for    each   interrogation; 
N  is  number  of   interrogation;    and 
1 2  is   twelve   unknown    for    4    transponders  coordinate. 

The    solution    of    Eguation    (A.  6)    is    N>3.      That   is,    the 
absolute   positions  in   three-dimensional   coordinates  of    all 
transponders   can    be    determined    with   as    few   as   three    consecu- 
tive   range    measurements  from   four   fixed   transponders. 

Mixhail   states  the  two-dimensional   distance    condition 
and    its    linearization    (Mikhail,    1981,    pp266-268),    which   can 
be    extented   to  three-dimensional  coordinates-    -  Equation 
(A.  5)    can  be   linearized; 

t=t°+l  (X.-U)  /Up  >X+[  (Y.-V)/iO]AY+[  (2-W)  /10]AZ-  '  (A. 7) 

I  C  L  Lit  (.(.(,  i  L 

The  adjusted  distance  is  then: 

5.       =     1-       +     V;  (A.  8) 

(.  t  u 

where        1-    are   observed   valuesof    the   distance;    and 
v-     are   corresponding   residuals. 

Eguation    (A. 7)    then    becomes; 

v-[   (X.-U)  /$o  ]AX.-[  (Y.-V)  /1°  ]AY.-[  (Z.-W)  /2?  >Z  =  io-J2.  (A  .  9) 

CO  £,  i-  L  LLC  CCC1 

Letting 

fL     =   JB.o    -    j?L  (A.  10) 

in    matrix   form    Equation     (A.  9)     becomes: 
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(A.  11) 


which  may  be  concisely  written  as: 


V     +  B    ax     =  f  t 


(A.  12) 


where 

V  is  residual  of  the  observation; 
B  is  numerical  coefficient   of  the  unknown 
variable  correction  ; 
A  x  is  unknown  variable  correction  ;  and 
f  is  numerical  constant  term  . 

The  weight  matrix  can  be  derived  from  the  standard  devi- 
ation of  the  slant  range,  which  discussed  in  the  previous 
chapter.   The  least  squares  solution  of  Equation  (A-  12)  is 
(Mikhail,  1981,  p259)  : 


N  =  B  KB 


(A.  13) 


t  =  B  Wf 


(A-  14) 


AX  =  N-i  t 


(A.  15) 


X  =  X0  +AX 


(A.  16) 


v  =  f  -  BAx 


(A.  17) 
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d-  =[2"2  (B  N"1  Bt)]»/2  (A.  18) 


2*2  =  v*  W  V  /  r  (A.  19) 


r  =  N  -  3  (A.  20) 

where 

N  is  coefficient  matrix  of  the  normal  equations, 
W  is  weight  matrix, 

t  is  vector  of  'constants'  in  the  normal  equations, 
Ax  is  vector  of  parameter  corrections, 
X0is  vector  of  approximate  value, 
X  is  vector  of  adjusted  value, 
d^is  standard  deviation  of  adjusted  distance, 
0-2  is  estimate  of  the  reference  variance, 
r  is  number  of  statical  degrees  of  freedom, 
N  is  number  of  given  transducer  position,  and 
3  is  number  of  observations  necessary  to  specify 

uniquely  the  model  that  underlies  the  adjustment 
problem. 

Since  during  linearization  we  neglect  all  second  and 
higher  order  terms,  we  must  ensure  that  X  is  not  signifi- 
cantly in  error  because  of  this  approximation. 
Consequently,  we  must  iterate  the  solution  by  using  X  as  a 
new  approximation  and  compute  anotner  correction.   This 
procedure  is  repeated  until  the  correction  is  insignificant. 
The  final  result  is: 


X„  =  x\  +  Ax  =  XM  (A.  21) 


wnere 
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n  is  total  iteration, 
£1  is  adjusted  value  at  n  iteration, 
Ax   is  correction  at  n  iteration  can  be  neglected  compare 
with  previous  correction  or  within  tolerance 
correction,  and 
X^  is  true  value  of  each  variable. 

Since  errors  in  transducer's  position   exist,  it  now 
must  be  established  how  this  error  affects  the  accuracy  of 
transponder  positions.   This  is  discussed  below: 

Let  the  standard  deviations  of  transducer  coordinates  be 

±6   ,  ±6,.    ±6  ,.   Kecall  Equation  (A.  5)  : 
u     v     W 

$.    =  [  (U  -X-  )2*  (V  -Y.  )  2  +  (»  -z.  )  2  ]i/2  (A.  22) 

^  C  (■  L» 

The    unknown  coordina tes. X, Y#Z    ,    are; 

X;   =     U    -    [J!?-  (V-  Y-  )   2-(W_Z.  )  2]l/2  (A.  23) 

Y-  =    V    -   [  i2-  (W-Z.  )  2-  (U-X.  )  2  ]i/2  (A.  24) 

Z.  =    W    -    [i#2-  (J-X.  )  2-(v-Y.  )  2]i/2  (A.  25) 
Now    assume 

F1    =   [12-(V-Y.  )  2-(W-Z.)  2]i/2  (A.26) 

F2    =   [$2-(W-Z.  )  2-(u-x.  )  2]i/2  <A.27) 

L  L     •      L 

F3  =  [f2-(u-x.  )  2-(V-Y.  )  2]i/2  (A. 28) 

So,  the  variances  of  Equations  (A.23),  (A. 24)  and  (A. 25)  are 
(Mikhail,  1981,  p. 181): 

fa     =d2-(i.  /F1)  2(<2«.  (v-YL/F  1)2(^2  «-  (W-Zj^/F  1)2  (J2  (A.  29) 
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6*    =C^2-(j?. /F2)  2(i2+  (W-ZL/F2)2dj+(U-X.u/F2)2d2  (A.30) 


6*    =d2-(i?. /F3)  2d2*  (U-X./F3)  2(^2+ (v-Y.  /F3)  2d2  (A. 31) 


and   the    position  accuracy   is: 


fT    =      ±    (d«    *dy2    +d2   )  1/2 


(A.  32) 
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